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Abstract 

Observations of comets and asteroids show that the Solar Nebula that spawned our 
planetary system was rich in water and organic molecules. Bombardment brought these 
organics to the young Earth’s surface, seeding its early chemistry [T]. Unlike asteroids, 
comets preserve a nearly pristine record of the Solar Nebula composition. The pres¬ 
ence of cyanides in comets, including 0.01% of methyl cyanide (CH3CN) with respect 
to water, is of special interest because of the importance of C-N bonds for abiotic amino 
acid synthesis [2]. Comet-like compositions of simple and complex volatiles are found 
in protostars, and can be readily explained by a combination of gas-phase chemistry 
to form e.g. HCN and an active ice-phase chemistry on grain surfaces that advances 
complexity [3]. Simple volatiles, including water and HCN, have been detected previ¬ 
ously in Solar Nebula analogues - protoplanetary disks around young stars - indicating 
that they survive disk formation or are reformed in situ^ O [HI [Z] • It has been hitherto 
unclear whether the same holds for more complex organic molecules outside of the 
Solar Nebula, since recent observations show a dramatic change in the chemistry at 
the boundary between nascent envelopes and young disks due to accretion shocks [8]. 

Here we report the detection of CH3CN (and HCN and HC3N) in the protoplanetary 
disk around the young star MWC 480. We find abundance ratios of these N-bearing 
organics in the gas-phase similar to comets, which suggests an even higher relative 
abundance of complex cyanides in the disk ice. This implies that complex organics ac¬ 
company simpler volatiles in protoplanetary disks, and that the rich organic chemistry 
of the Solar Nebula was not unique. 
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MWC 480 is a Herbig Ae star with an estimated stellar mass of 1.8 Solar masses (Mq) 
[2] in the Tanrus star forming region at a distance of 140 pc. The star is surrounded 
by a 0.18±0.1 Mq protoplanetary disk, an order of magnitude more massive disk than 
the Minimum Mass Solar Nebula[T 0 l |TT]. Compared to disks around Solar-type stars, the 
MWC 480 disk is 2~3 times warmer at a given radius [T^ 113] and exposed to orders of 
magnitude higher ultraviolet radiation levels. Despite these environmental differences the 
volatile composition in the MWC 480 disk appears largely similar to that of disks around 
Solar-type stars, except for lower abundance of cold (T<20 K) chemistry tracers in the 
outer disk m I15j . The inner disk chemistry of MWC 480 has not been studied, but 
Herbig Ae and T Tauri disks are observed to present different volatile compositions close 
to their stars [El- 

Using the Atacama Large Millimeter/submillimeter Array (ALMA), we detected two 
emission lines of CH 3 CN from the MWC 480 protoplanetary disk (the 14o — 13o line at 5a, 
and the 14i — 13i line at 3.Sc). We also detected emission lines from the N-bearing carbon 
chain HC 3 N, and from the 13-C isotopologue of HCN. We targeted the 13-C isotopologue 
rather than the more abundant H^^CN because the latter is optically thick and therefore 
a poor tracer of the HCN abundance. Figure 1 shows the spatially resolved line detec¬ 
tions together with a dust continuum emission map, demonstrating the spatial coincidence 
between CH 3 CN, HC 3 N, HCN and dust emission from the disk. The angular resolution 
is 0.”4 — 0”.6, corresponding to 50-70 AU (1 AU is the distance from the Earth to the 
Sun). The emission is also spectrally resolved, which can be used to probe smaller spatial 
scales. Figure 1 shows the velocity gradient across the disk that arises from Keplerian 
rotation in all three lines. Figure 2 shows the spectra of the three lines; each displays the 
double-peaked structure typical of a rotating disk. Table 1 lists the integrated line fluxes. 

We use the spectrally and spatially resolved line emission to constrain the radial profiles 
of molecular column density based on parametric abundance models defined with respect 
to the adopted MWC 480 disk density and temperature structure (Methods). Figures 2 
and 3a-c show the synthetic line spectra and maps from the best-fit models, demonstrating 
the good match between the models and data. Figures 3d-f shows the best-fit radial column 
density profile, together with all profiles consistent with the data within 3a, and abundances 
at 30 and 100 AU are reported in Table 1. The best-fit profiles have different slopes for the 
different molecules. The H^^CN column density decreases with radius, which is consistent 
with predictions from disk chemistry models [IT] . The increasing column density with 
radius of HC 3 N, effectively a ring, is not predicted by models [18]. indicative of that disk 
chemistry models are incomplete for HC 3 N. The CH 3 CN emission is best reproduced with 
a flat profile, but other profiles cannot be excluded. 

The absolute abundances depend on the assumed disk density structure, and we there¬ 
fore compare the complex cyanides in the MWC 480 disk to comet and to protostellar com¬ 
positions using HCN as a reference species, similar to the practice in cometary studies [19]. 
We calculate abundances of CH 3 CN and HC 3 N with respect to HCN at 30 AU, the smallest 
disk radius accessible by the ALMA observations, and at 100 AU, the outer boundary of 
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the cyanide emission maps. Acconnting for the higher luminosity of MWC 480 compared 
to the young Sun, this radial range in the MWC 480 disk corresponds to the comet forming 
zone of 10-30 AU in the Solar Nebula[T9j. Assuming a standard HCN/H^^CN ratio of 70, 
the best-fit HC 3 N and CH 3 CN abundances with respect to HCN are 0.4 and 0.05 at 30 AU, 
and 5 and 0.2 at 100 AU, respectively. The CH 3 CN/HCN abundance ratios are robust to 
model assumptions to within factors of a few, while the HC 3 N/HCN abundance ratio may 
be overestimated by an order of magnitude using our simple abundance model (Methods). 
Conservatively, both the CH 3 CN and HC 3 N abundances with respect to HCN are thus 
~5% at 30 AU and ~20% at 100 AU. A typical comet contains 10% of CH 3 CN and HC 3 N 
with respect to HCN[T9]. The MWC 480 gas-phase cyanide composition at both 30 and 
100 AU is thus cometary within the observational and model uncertainties. 

The relationship between gas-phase abundance ratios and the abundance ratios in ices, 
the main volatile reservoirs in disks pUl 1171118] . depends on both desorption characteristics 
and chemistry (Methods). HCN, HC 3 N and CH 3 CN are characterized by similar freeze-out 
and desorption kinetics, but different chemical pathways. In particular, the existence of 
efficient grain surface formation pathways to CH 3 CN enhances CH 3 CN with respect to 
the other cyanides in the ice mantles. The scale of this enhancement factor varies among 
models, but it is at least one order of magnitude (Methods). This results in an expected 
minimum CH 3 CN/HCN ice ratio of 0.5 at 30 AU in the MWC 480 disk, considerably higher 
than the typical Solar System comet. 

The CH 3 CN/HCN and HC 3 N/HCN ratios in the MWC 480 disk are also high when 
compared to protostars. The HCN/HC 3 N/CH 3 CN ratio is 1/0.01/0.08 toward the Solar- 
type protostellar binary IRAS 16298-2422|21j. and similar abundance ratios are found to¬ 
ward more massive systems. The MWC 480 cyanide composition is thus difficult to explain 
by inheritance alone, as has been suggested for e.g. H 20 [ 7 ]- Rather the observed high disk 
abundances likely reflect an efficient disk chemistry that readily converts a large portion of 
the carbon originally in CO and other small molecules into more complex organics |22 1 123j 
during the first million years of the disk life time. 

This early, efficient complex chemistry in protoplanetary disks impacts the surface con¬ 
ditions of rocky (exo-)planets. In the Grand Tack model of early Solar System dynamics |24j. 
the coupled migration of Jupiter and Saturn in the gas-rich disk from which they formed 
caused scattering of volatile-rich (and thus complex organic-rich) planetesimals inwards, 
mixing with the asteroid belt, and outwards, producing the present-day Kuiper belt and 
comets [25]. While comets are expected (and observed) to conserve most of their original 
compositions due to low temperatures and thus very long desorption and chemical time 
scales [T^. asteroid belt bodies are more processed. A combination of evaporation and 
chemistry probably resulted in both a net loss of volatiles and in a relative enhancement 
of complex organics with respect to H 2 O over the subsequent 10s of Myr[26|. A second 
instability 100 million years later resulted in a heavy bombardment of the Earth by mainly 
asteroid belt bodies, including the icy bodies originating in the outer Solar Svstem[27|. 
The frequency of these instabilities in exo-planetary systems is unknown, but the pile-up 
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of Hot Jupiters and Super-Earths reveal that major planetary migrations are common. 
When these migrations cause bombardment of icy bodies, they are likely organic-rich. 

The high ratio of complex to simple cyanides in the MWC 480 disk implies that the rich 
organic composition of comets is not unique to our Solar System, and could be common¬ 
place. Laboratory experiments have shown that the same ice chemistry that produce 
CH 3 CN, i.e. photo processing of interstellar ice analogs also produce simple sugars and 
amino acids [28l I29j . This suggests that the early surface conditions of Earth, set by comet 
and asteroid bombardment, may be common for young rocky planets, and that conditions 
favourable to an even richer chemistry may be ubiquitous. 
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Table 1: Molecular data. 


Molecule 

QN 

Rest freq. 

Eu 

Integrated 
flux [rms] 

•^lOOAU 



GHz 

K 

mJy km/s ( 10 “^^ 

(10-13 n^i) 


beam ^ 


H13CN 

J=3-2 

259.0118 

24.9 

47[8] 

1.6 [1.3-2.0] 

1.6 [1.3-2.0] 

HC 3 N 

J=27-26 

245.6063 
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66 [ 8 ] 

43 [36-54] 

480 [400-600] 

CH 3 CN 

14o - 13o 

257.5274 

92.7 

30 [ 6 ] 

4.8 [2.3-7.5] 

16 [6-55] 

CHgCN 

14i - 13i 

257.5224 

100 

23 [ 6 ] 




^best-fit abundances assuming a vertically constant abundance, the 3a abundance range 

is in []. 
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Figure 1: ALMA detections of simple and complex cyanides in the MWC 480 
protoplanetary disk, a, 1.1 mm emission (black contours are 3cr + a x b-d, 

integrated emission of H^^CN (a), HC 3 N (b) and CH 3 CN (c) lines (colour: see colour scale 
on the right). Black contours are [3,4,5,7,10]cr. e-g, as b-d, but for 2 km/s velocity bins 
around the source mean velocity, displaying the disk rotation. Positions are relative to 
the continuum phase centre (marked with a cross) at right ascension (RA) 04:58:45.94 and 
declination (DEC) +29:50:38.4. The synthesized beam is shown in the bottom left corner 
of each panel. 
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Figure 2: Spectra of detected cyanides in the MWC 480 protoplanetary disk. 

The observed spectra (contours) of H^^CN, HC 3 N, and CH 3 CN are extracted from ALMA 
spectral-image data cubes. The synthetic spectra (histograms) are based on the best fit 
disk abundance models in Fig. 3. The CH 3 CN spectrum contains two partially overlapping 
lines identified with the 14o — 13o and 14i — 13i transitions. The spectra were extracted 
from the spectral image cubes using a Keplerian mask to maximize the signal-to-noise. 
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Figure 3: Models of cyanide emission and radial distributions in the MWC 480 

disk, a-c, synthetic ALMA observations without added noise based on the best-fit radial 
column density profiles of H^^CN, HC 3 N, and CH 3 CN (colour: see colour scale on the 
right). Black contours are as in Fig. 1. d-e, best-fit radial column density profiles (thick 
dark blue line), shown together with all column density profiles that are consistent with 
data within 1, 2 and 3cr confidence intervals. Increasing deviations from the best-fit are 
shown with successively lighter shades and thinner lines. 
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Methods 


ALMA observations and data reduction. 

MWC 480 was observed with the Atacama Large Millimeter/submillimeter Array as a part 
of the Cycle 2 proposal 2013.1.00226. Observations were carried out on 2014 June 15, with 
33 antennas and baselines of 18-650 m (15-555 kA). The total on source integration time 
was 22 minutes. The nearby quasar J0510+1800 was used for phase and gain, bandpass and 
absolute flux calibration (J0510+1800 was itself calibrated to an absolute flux of 1.50[0.08] 
Jy on 2014 June 29). The correlator was configured to observe 14 spectral windows (SPWs) 
with resolution 6i'= 61.04 kHz. The H^^CN 3-2 line was placed in SPW 1, HC 3 N 27-26 in 
SPW 9 and CH 3 CN 14o-13o and 14i-13i in SPW 5, each 60-120 MHz wide and centered 
on 258.983, 245.622 and 257.498 GHz, respectively. 

The visibility data were calibrated by ALMA/ NAASC staff following standard proce¬ 
dures. Each SPW was further self-calibrated in phase and amplitude in CASA 4.2.2. The 
self-calibrated phase and gain SPW-specific solutions were applied to create individual 
SPW line data cubes using the Briggs weighting scheme, with Brigg’s robustness param¬ 
eter set to 1-2. For the weak CH 3 CN lines, a taper of 1" was applied to the visibilties 
to maximize the signal-to-noise. The data cubes were subsequently continuum subtracted 
using channels identified as free of spectral line emission. The resulting rms noise in 1 
km/s bins was ~2.3 mJy per beam. A CLEAN deconvolution was performed with a mask 
produced manually based on the H^^CN emission and applied to the HC 3 N and CH 3 CN 
spectral line cubes (after verifying that the HC 3 N and CH 3 CN cubes contained no addi¬ 
tional emission). For the 257 GHz continuum map, obtained by combining all SPWs in 
the USB), the synthesized beam is 0."65 x 0.''39 (PA=-7°) and rms noise 0.3 mJy beam“^. 
The 257 GHz continuum flux density was determined to be 331±33 mJy (assuming a 10% 
uncertainty in the absolute calibration), which agrees well with previous Submillimer Array 
results^^. 

Integrated emission maps were produced using the immoments task in CASA, adding 
emission >1(7 from each channel. Lower clipping levels do not change the result at the 
10% level. Spectra were extracted using both the empirical CLEAN masks and theoretical 
masks based on Keplerian rotation to isolate disk emission in each channel. The CLEAN 
masks resulted in the best signal-to-noise ratio, and these results are used in the study. 

Parametric disk density and temperature model. 

We adopt a previously developed disk density structure of MWC 480[T0]. The total surface 
density is given by the similarity solution, which self-consistently describes the density 
structure of an accretion disk^^’^^: 



( 1 ) 
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where Sc = Mgas{‘i — 7)/(27ri?c) is a normalization coefficient, Rc is the characteristic 
radius, and 7 is a gradient parameter describing the density fall-off with radius. The 
two-dimensional density structure in cylindrical coordinates is 

where H{r) = HiQQ{r/1Q{)AU)^ is the scale height, and h a parameter that describes the 
amount of flaring of the disk. The density structure parameters were constrained using 
high angular resolution observations of the MWC 480 dust continuum emission at 1.3 mm 
and 3 mm wavelengths m and are listed in Extended Data Table 1. The resulting mass 
surface densities in g cm“^ are 550, 91 and 4.7 at 1, 10 and 100 AU, respectively. 

Based on existing model fits to dust and observations^^’^^, the disk midplane 

temperature is parametrized as T{r) = Tq > where To is set to 23 K at = 100 AU, 

and the power-law index q is 0.4. To account for the presence of a vertical temperature 
gradient due to heating of the disk surface by the central star^^’^®, we parameterize the 
two-dimensional temperature structure in cylindrical coordinates as 

where Tq is the midplane temperature at Tq, g is a power-law index describing the de¬ 
crease in midplane temperature with radius, and /3 is a factor describing the increase in 
temperature at increasing disk height, which is set to 1.5^'^. The density and temperature 
disk structures are shown in Extended Data Fig. 1. 

Molecular abundance retrieval. 

To retrieve molecular abundance and column density profiles, we define parametric abun¬ 
dance models with respect to the adopted MWC 480 disk density and temperature struc¬ 
ture. We use a simple power-law prescription x = xiooAU x (r/100 AU)°‘ for the molecular 
abundances where x is the abundance of the molecule with respect to the total hydrogen 
density, xiooAU is the abundance at 100 AU, r is the disk radius in AU, and a is a power 
law index^®. To obtain a column density profile, which is more intuitively related to ob¬ 
servations, the abundance prescription is multiplied by disk surface density profile, Eq. 1. 
We set an outer cut-off radius of 100 AU, corresponding to the most extended emission 
observed for any of the cyanides. We also explored Tout of 60 and 80 AU, but this did 
not improve the fit for any of the molecules. For each molecule, we calculate a grid of 
abundance models that covers xiooAU from 10 “^^ — 10 “^° with respect to nn, and a be¬ 
tween 0 and 1 for H^^CN, 1 and 2 for HC 3 N, 0 and 2 for CH 3 CN. These initial ranges of a 
were selected based on visual inspection of the observed emission, noting its level of central 
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concentration, i.e. a = 0 corresponds to a steeply decreasing column density profile, a = 1 
to an almost flat profile, and a = 2 to an increasing profile with radius. 

The best-fit models are obtained by minimizing the weighted difference between the 
real and imaginary part of the complex visibility measured in the (u, u)-plane sampled by 
the ALMA observations^®’^'^’^®. We use the three-dimensional Monte Carlo code LIME^® 
to calculate the radiative transfer and molecular excitation. In addition to the molecular 
abundance profile and the disk density and temperature structure parameters, the radia¬ 
tive transfer modelling requires information on disk inclination, turbulence and velocity 
field. We adopted a disk inclination of 37°, disk turbulence of 0.05 km/s and a Keplerian 
velocity field with a stellar mass of 1.8 Mq m- For HCN and HC 3 N, the level populations 
were computed in non-LTE (Local Thermodynamic Equilibrium) with collision rates listed 
in the BASECOL database^^’^^. For CH 3 CN, we assumed LTE, which is a reasonable ap¬ 
proximation, since the expected disk emissive layers have higher densities than the critical 
density. We checked this assumption by running several HCN and HC 3 N models assuming 
LTE and found a good agreement (within 20%) between LTE and non-LTE models. In 
Extended Data Fig. 2 we show the best-fit models for different a. In the case of HC 3 N, 
only the models with a = 2 shows the observed ring-like structure and are therefore the 
only ones considered further. 

Within the adopted model framework, the molecular abundances and abundance ratios 
are constrained within a factor of a few. These models implicitly assume a constant abun¬ 
dance profile with disk height, which may be a poor approximation in the outer disk where 
freeze-out in the midplane results in a vertically layered structure[2n|’^^’^^. The existing 
data are, however, insufficient to put any constraints on the vertical distribution. Rather 
than pursuing a more complex model, we opt to test the sensitivity of the derived abun¬ 
dances and abundance ratios on the details of the vertical structure assumptions. We ran 
a second grid of models to simulate the effect of freeze-out in the midplane by reducing the 
molecular abundances by three orders of magnitude at zjR < 0.2. Using the same fitting 
procedure as before, we find best-fit xiooau of 1 x 10“^^ (a=0), 3.5 x 10“^*^ (a=2), and 
5 X 10“^^ (a=l), for H^^CN, HC 3 N and CH 3 CN, respectively. These abundances are an 
order of magnitude higher than those obtained assuming a vertically constant abundance. 
The CHsCN/H^^CN abundance ratio profile is unchanged within the model uncertainties, 
however. Despite a higher excitation level of CH 3 CN the two molecules emit from a sim¬ 
ilar disk layer when assigned the same vertical abundance profile due to a combination of 
a high critical density of the CH 3 CN transition and a rapid fall-off in density with disk 
height, i.e. practically no CH 3 CN emission originates at zjR > 0.5. In contrast, the de¬ 
rived HC 3 N/HCN ratio does depend on the vertical structure. The HC 3 N transition has 
an even higher energy level and a lower critical density than CH 3 CN, which implies that 
some of the HC 3 N emission can originate in more elevated disk layers than H^^CN and 
CH 3 CN. The reported CHsCN/H^^CN results are thus robust to the abundance model 
assumptions, while the reported HC 3 N/H^^CN results may be overestimated by up to an 
order magnitude. 
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Disk chemistry modelling: constraining the ice-to-gas ratios. 

To constrain the origin of CH 3 CN and therefore how its ice-to-gas ratio compares with 
HCN, we first explore how much of the observed CH 3 CN gas could come from gas-phase 
chemistry. We ran a grid of pseudo-time dependent models with a complete gas-phase 
chemistry*^^, but without any grain surface reactions (except for adsorption and desorp¬ 
tion, and H 2 formation on grains), for 1 Myr. The grid covers several orders of magnitude 
of densities, temperatures, UV fields and ionization fractions - the four most important 
regulators of gas-phase chemistry in disks. For gas-phase CH 3 CN production, the ion¬ 
ization rate is expected to be of special importance because the main gas-phase forma¬ 
tion pathway includes CH(|' (through HCN-|-CH^, which has a expert-validated reaction 
rate in the KID A database^^). Extended Data Figure 3 shows that none of these models 
can produce CH 3 CN/HCN >0.01 and that in the vast majority of the parameter space 
CH 3 CN/HCN< 0 . 001 . a ratio of 0.01 is approached for an ionization rate of 10“^^ s“^, but 
such high ionization rates are only attained in the disk atmosphere where stellar X-rays 
regulate the ionization balance. This disk layer contain a very small fraction of the total 
disk mass and therefore contributes a negligible amount of the total molecular column. 
Deeper into the disk, the ionization rate is reduced below 10“^^ s“^ because of attenuation 
of both X-rays and cosmic rays^®. A ratio of 0.01 is also approached at lower ionization 
rates at very low temperatures and high densities. Such environments are characteristic for 
the outer disk midplane, which contain a lot of mass, but freeze-out at these temperatures 
results in CH 3 CN abundances below 10“^^, which would not contribute to the observed 
CH 3 CN emission. The observed high abundance of CH 3 CN gas with respect to HCN there¬ 
fore implies that grain surface chemistry contributes significantly to the observed CH 3 CN 
abundance. 

We ran a number of complete disk chemistry models for IMyr, the estimated age 
of the MWC 480 star-|-disk system, which include gas and grain surface chemistry, as 
well as different levels of turbulence, to calculate the range of plausible ice-to-gas conver¬ 
sion factors for HCN (and other gas-phase chemistry products) and CH 3 CN^'^. The main 
grain surface CH 3 CN formation pathways in these models are hydrogenation of C 2 N and 
CH 3 -I-CN and ice photochemistry. The main desorption pathway is UV photodesorption 
because of the high binding energies of HCN. HC 3 N and CH 3 CN: 4170, 4580, and 4680 K, 
respectivelyrespectively, corresponding to sublimation temperatures of ~90 K at a 
gas density of 10® cm“^ and 110 K at the gas density of 10^® cm“^. The three cyanides 
should also present very similar UV photodesorption efficiencies because of similar radia¬ 
tion cross-sections and binding energies®^. Absent grain surface chemistry, these molecules 
would thus have had comparable ice-to-gas ratios. 

Because CH 3 CN can form on the grains in the cold midplane a key feature is the 
coupling between the ice formation layer and the disk layers where molecules are efficiently 
desorbed into the gas-phase. Without such a coupling, ice chemistry cannot affect the gas- 
phase volatile composition. In the disk models, we define the vertical diffusion coefficient 
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as Dz = azC^I^, where is a free parameter, Cg is the local sound velocity, and Q is 
the Keplerian orbital frequency. Different levels of turbulence are simulated by different 
values of Oz between 0 and 0.01. In models with turbulence turned off (a^ = 0), there 
is no mixing between the cold, icy midplane and the intermediate layers where gas-phase 
molecules are abundant (and the observed cyanides likely reside). Extended Data Figure 
4a shows that without turbulence the predicted CH 3 CN abundances are low everywhere 
(< 10“^^ w.r.t tt-h, where uh is the number density of H nuclei) and the vertically averaged 
CH 3 CN/HCN abundance ratio is <0.001; see Extended Data Figure 5a,d). 

In models that include vertical mixing (a^ = 0.001 — 0.01), the CH 3 CN abundance as 
well as the CH 3 CN/HCN abundance ratio are enhanced by an order of magnitude because 
of mixing of midplane icy grains into UV exposed disk regions where the ices photodes- 
orb (Extended Data Figs. 4 and 5). In particular the CH 3 CN/HCN ratio approaches 
10“^ outside of 30 AU, within an order of magnitude of the observed ratios. Model-data 
comparison thus suggest that the intermediate disk layers are strongly coupled to the icy 
midplane through mixing. This implies that gas-phase observations of CH 3 CN and other 
cyanides at intermediate disk heights can be used to trace the total volatile reservoir, even 
if the conversion from gas to ice can be complex. 

Using the three models, with no, moderate and strong mixing we calculate the HCN and 
CH 3 CN ice-to-gas density and column density ratios as a function of disk radius (Extended 
Data Fig. 6,7). For HCN, the ice-to-gas ratio varies between 10® and 10® between 30 and 
100 AU. For CH 3 CN the ice-to-gas ratio varies between 10^ and 10“^. In other disk chemistry 
models [T 8 ]’®^ the contrast between the two ice-to-gas ratios is even higher. Considering the 
range of values and their sensitive dependence on assumptions about disk turbulence, we 
conclude that the ice-to-gas ratio for CH 3 CN is at least 10 times higher than HCN and 
HC 3 N. At longer time scales the balance between formation in the midplane and destruc¬ 
tion at more exposed disk layers determines whether the CH 3 CN abundance increases or 
decreases - in general the more turbulence the faster the destruction. 

It is important to note, that we employ generic disk structures®® that have not been ht to 
the MWC 480 data. The temperature and density structures are comparable (cf. Extended 
Data Figs. 1 and 4). The UV held is lower than would be expected toward a Herbig Ae 
star and the gas-phase CH 3 CN production may therefore be somewhat overestimated in 
our model, but as shown above, gas-phase chemistry is a minor contributor to the overall 
CH 3 CN gas budget. The grain surface production takes place in UV shielded regions and 
should thus not be affected by an increased UV flux. An increased UV held intensity 
may change the desorption/adsorption balance, however, since UV photodesorption is an 
important desorption mechanism. This increase in photodesorption will be accompanied 
by an increased level photodissociation and the main effect will therefore be to push the 
desorbing layer closer to the midplane rather than increasing its abundance. An A star 
type radiation held may aid the detection of gas-phase CH 3 CN, since less turbulence may 
be required to bring the icy grains up to UV exposure and into the gas-phase. A decrease 
in the disk height of the cyanide gas layer, compared to our generic disk model, may 
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also reduce the absolute ice-to-gas ratios, but the relative ratios of species with similar 
adsorption and desorption characteristics should not be significantly affected. 

Desorption Mechanism 

In our disk chemistry model, UV photodesorption regulates the desorption of oh HCN, 
HC 3 N and CH 3 CN ice into the gas-phase. This result depends on a combination of disk 
and star properties and on the theoretical yields of different desorption mechanisms, many 
of which are poorly constrained. While relative ice-to-gas ratios of related species are 
insensitive to the details of the desorption process, the absolute ice-to-gas ratios found in 
our models must be considered highly uncertain. Based on theory, different desorption 
mechanisms should present different radial and vertical abundance distributions, which 
could be tested observationally. Desorption due to release of chemical formation energy 
should be most important in the cold disk midplane, while UV photodesorption is only 
efficient in warmer, upper disk layers, and thermal desorption in the disk layers that are 
warm enough. In particular, colder molecular emission layers should be characterized by 
lower excitation temperatures, which can be measured when multiple lines of the same 
species are observed. Observational constraints on ice desorption in disks will thus be 
attainable with ALMA. 

Different desorption mechanisms should also exhibit different chemical signatures, al¬ 
though the details are often model-dependent. In the case of thermal desorption, the gas 
phase composition in a layer should simply depend on volatility, albeit with the added 
complication of ice entrapment^^. Chemistry can also help distinguish between ice and 
gas origins and thus provide an independent measure of the relative importance of the 
two pathways, he high ratio of CH 3 CN/CH 3 NC in the Horsehead PDR (Photon Domi¬ 
nated Region), for example, was used as evidence for a grain surface origin of the observed 
CH 3 CN®^. Chemical tracers could also provide useful constraints on the incident radiation 
fields and temperature structure and thus aid in ruling out some of the potentially impor¬ 
tant desorption mechanisms. For example CN/HCN is a proposed tracer of UV flux and 
HNC/HCN of temperature^®’®'^. 
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Extended Data 


Extended Data Figure 1: Physical model for the disk of MWC 480 


Parameters 

Values 

Stellar properties 


Estimated distance: d (pc) 

140 

Stellar mass: M* (Mq) 

1.8 

Disk structure properties 

Disk mass: M(i(MQ) 

0.18 

Characteristic radius: Rc(AU) 

81 

Outer cut-off radius (AU) 

100 

Scale height: -ffiooAU (AU) 

16 

Flaring index: h 

1.25 

Density power-law index: 7 

0.75 

Midplane temperature: Tiooau (K) 

23 

Temperature power-law index: q 

0.5 

Vertical temperature gradient index: /3 

1.5 

Disk geometric and kinematic properties 

Inclination: i (deg) 

37 

Systemic velocity: Vi;, 5 j?(kms“^) 

5.0 

Turbulent line width: Vj„rf,(km s“^) 

0.05 

Position angle: P.A.(deg) 

58 
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Extended Data Figure 1: Model of the physical structure of the MWC 480 proto¬ 
planetary disk, a, radial and vertical disk temperature profile (color: see color scale on 
the right, contours: 20,30,50,100, and 1000 K) with important temperatures marked with 
contours, b, radial and vertical density profile (color: see color scale on the right, contours: 
10^*^, 10®, 10® and 10^ cm“®). The disk height over radius of 0.2 is marked as well. 
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Extended Data Figure 2: Synthetic observations of H^^CN, HC 3 N and CH 3 CN. 

The models are based on best-fit to data for different choices of a. a-g, integrated emission 
maps (colour: see colour scale on the right). Black contours are the observed [3,4,5,7,10]cr 
in Fig. 1. The synthesized beam is shown in the bottom left corner of each panel. Note 
the change in emission profile between a = 1 and 2 for HC 3 N. 


20 













UV radiation field 





100 150 200‘ 



50 100 150 200 ‘ 

Temperature [K] r- 


50 100 150 200 

Temperature [K] 


10 '” 

10 * 

10 ° 

10 ^ 



\ 

jLn 

O' 

e 



10 '° 

10 ° 

10 ° 

10 ° 




f 



Temperature [K] 


50 100 150 200 ‘ 

Temperature [K] 


50 100 150 200 

Temperature [K] 





50 100 150 200‘ 

Temperature [K] 


50 100 150 200 ‘ 

Temperature [K] 


50 100 150 200 

Temperature [K] 


10 '“ 

10 ° 

10 ° 

10 ° 


V /* i -'v 

'g 10'° 
JL 10 ° 
^ 10 ° 

■./ 

'g 10 ’° 

JL 10 ° 
^ 10 ° 

J y" : 

j_^^^_ 

2 10 ° 

k , _ 

1 10 ° 

1 .—-5- 


50 100 150 200 

Temperature [K] 


50 100 150 200 

Temperature [K] 


50 100 150 200 

Temperature [K] 


-5 -4 -3 

log(CH 3 CN/HCN) 


Extended Data Figure 3: Models of gaseous CH 3 CN/HCN abundance ratios un¬ 
der different physical conditions, a-1, the CH 3 CN/HCN abundance ratio on a loga¬ 
rithmic scale (colour: see colour scale on the bottom). UV radiation flux increases from 
left to right from Gq = 1 (a,d,g,j) to Gq = 10 (b,e,h,k) to Gq = 100 (c,f,i,l)) where Gq is 
scaling factor in multiples of the local interstellar radiation field. The ionization rate of H 2 
increases from top to bottom from 10 “^^ s“^ to 10 “^® s“^ (d^f) to 10 “^^ s“^ (g^i) 

to 10 “^^ s“^ 
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Extended Data Figure 4: Models of gaseous CH 3 CN in disks with and withont 
turbnlent diffusion, a, the abundance of CH 3 CN (colour: see colour scale on the right) 
as function of disk radius {R) and height scaled by the radius {ZjR) in a model without 
turbulence. The dashed lines indicate gas temperatures of [30,50,100] K. b-c as a, but in 
disk models that include turbulence parameterized by = 10 “^ (b) and = 10 “^ (c). 
d, the vertically integrated column density of CH 3 CN from a-c (solid line: = 0, dashed 
line: = 10 “^, dotted line: = 10 “^). 
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Extended Data Figure 5: Models of gaseous CH 3 CN/HCN ratios in disks with 
and without turbulent diffusion, a-d, as in Extended Data Fig. 4. 
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Extended Data Figure 6 : Models of gas-to-ice ratios of HCN in disks with and 
without turbulent diffusion, a-d, as in Extended Data Fig. 4. 


24 


log(lce/gas abundance ratio) 














I 


5 


4 


3 


2 


Extended Data Figure 7: Models of gas-to-ice ratios of CH 3 CN in disks with and 
withont tnrbulent diffusion, a-d, as in Extended Data Fig. 4. 
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